First Principles Studies of Multiferroic Materials 

Silvia Picozzi 1 and Claude Ederer 2 

1 Consiglio Nazionale delle Ricerche - Istituto Nazionale per la Fisica della Materia 
(CNR-INFM), CASTI Regional Laboratory, 67100 L'Aquila, Italy 
2 School of Physics, Trinity College, Dublin 2, Ireland 

Abstract 

Multifcrroics, materials where spontaneous long-range magnetic and dipolar orders coex- 
ist, represent an attractive class of compounds, which combine rich and fascinating funda- 
mental physics with a technologically appealing potential for applications in the general area 
of spintronics. Ab-initio calculations have significantly contributed to recent progress in this 
area, by elucidating different mechanisms for multifcrroicity and providing essential informa- 
tion on various compounds where these effects are manifestly at play. In particular, here we 
present examples of density-functional theory investigations for two main classes of materials: 
a) proper multiferroics (where ferroelectricity is driven by hybridization or purely structural 
effects), with BiFe03 as prototype material, and b) improper multiferroics (where ferroelec- 
tricity is driven by correlation effects and is strongly linked to electronic degrees of freedom 
such as spin, charge, or orbital ordering), with rare-earth manganites as prototypes. As for 
proper multifcrroics, first-principles calculations are shown to provide an accurate qualitative 
and quantitative description of the physics in BiFe03, ranging from the prediction of large 
ferroelectric polarization and weak ferromagnetism, over the effect of epitaxial strain, to the 
identification of possible scenarios for coupling between ferroelectric and magnetic order. 
For the class of improper multiferroics, ab-initio calculations have shown that, in those cases 
where spin-ordering breaks inversion symmetry (i.e. in antiferromagnetic E-type H0M11O3), 
the magnetically-induced ferroelectric polarization can be as large as a few /zC/cm 2 . The 
presented examples point the way to several possible avenues for future research: On the 
technological side, first-principles simulations can contribute to a rational materials design, 
aimed at identifying spintronic materials that exhibit ferromagnetism and ferroelectricity at 
or above room-temperature. On the fundamental side, ab-initio approaches can be used to 
explore new mechanisms for ferroelectricity by exploiting electronic correlations that are at 
play in transition metal oxides, and by suggesting ways to maximize the strength of these 
effects as well as the corresponding ordering temperatures. 



1 Introduction to multiferroic materials 

Recent years have seen an enormous increase in research activity in the field of multiferroic 
materials and magneto-electric effects. In December 2007 Science Magazine listed multiferroic 
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materials as one out of ten "Areas to watch in 2008", the only entry from the Materials Sci- 
ence/Condensed Matter area that was included in this list. First principles calculations using 
density functional theory (DFT) [H HI [3] have played an important role in this "Renaissance of 
Magnetoelectric Multiferroics" [1] . In the present paper we give a brief summary of the current 
status of research on multiferroic materials and highlight some of the contributions that have 
been made using first principles electronic structure calculations. 

According to the original definition put forward by Hans Schmid |5j, multiferroic materials are 
materials that combine two or more of the primary forms of ferroic order, i.e. ferroeleasticity, 
ferroelectricity, ferromagnetism, and ferrotoroidicity. In practice, most of the recent research has 
focused on materials that combine some form of magnetic order (ferromagnetic, antiferromag- 
netic, non-collinear, . . . ) with ferroelectricity. Therefore, the term multiferroics is nowadays 
often used synonymous with magnetic ferroelectrics. 

Research on multiferroics (or magnetic ferroelectrics) is also intimately interwoven with research 
on the magneto- electric effect, which is the property that in certain materials a magnetic field 
induces an electric polarization and, conversely, an electric field induces a magnetization. Tra- 
ditionally, one distinguishes between linear, quadratic, and higher order magneto-electric effects 
[B], but more recently the term "magneto-electric effect" is often (mis-)used to describe any 
form of cross-correlation between magnetic and (di-)electric properties. (For example, when the 
application of an external magnetic field induces a phase transition between ferroelectric/non- 
ferroelectric phases.) It is important to point out, though, that not every magnetic ferroelectric 
exhibits a linear magneto-electric effect (in the original sense) and that not every material that 
exhibits a linear magneto-electric effect is also simultaneously multiferroic. 

Due to the combination of magnetic and dielectric properties, with eventual cross-coupling be- 
tween these properties, multiferroics have immense potential for technological device applications 
and at the same time they pose very interesting and rich fundamental physics problems. It is 
probably this combination of applied and fundamental research that is partly responsible for the 
strong attraction that these materials have developed in recent years. 

Multiferroics form a very diverse class of materials, and there is no unique "theory of multi- 
ferroics". Nearly every material has to be studied on its own right, and eventually involves 
very different physical mechanisms than other multiferroic materials. However, it has proven 
to be very useful to classify different multiferroics according to the mechanism that drives the 
ferroelectricity in the corresponding systems. In particular two major classes of multiferroics 
can be distinguished: 

1. Multiferroics, where the ferroelectricity is driven by hybridization and covalency or other 
purely structural effects. 

2. Multiferroics, where the ferroelectricity is driven by some other electronic mechanism, e.g. 
"correlation" effects. 

In the second case, ferroelectricity always arises as a secondary effect that is coupled to some 
other form of ordering, such as magnetic or charge ordering. Therefore, these systems are 
often called "improper magnetic ferroelectrics". We note that also in the first class at least 
one material, hexagonal YMnC>3, has been classified as an improper ferroelectric, where the 
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electric polarization is not the primary order parameter, but instead is coupled to a different 
non-polar structural instability [7j. In spite of that, and for the purpose of this article, we 
will call materials belonging to the first category "proper magnetic ferroelectrics" (or "proper 
multiferroics" ) , whereas materials in the second category will be called "improper magnetic 
ferroelectrics" (or "improper multiferroics" ) . We note that this "zoology" of multiferroics is still 
work in progress, and that the discovery of new materials might require a further refinement or 
redefinition of previous classifications. 

In this article, we are not attempting to provide a complete review of all first principles work 
that has been carried out so far. Instead, we discuss some specific examples that illustrate the 
power of these methods in elucidating the physical origins of the observed properties of known 
multiferroics, and point out the possibilities in predicting novel effects and designing new mate- 
rials with optimized properties. Also, we focus only on single-phase (bulk) materials, therefore 
leaving out all those effects coming from the combination of ferroelectrics and ferromagnets in 
(artificial) multiferroic heterostructures. Several excellent review articles about general aspects 
of multiferroic materials and magneto-electric effects have already been published, see for exam- 
ple Refs. (H llHl El [121 [13] , and much of the early first principles work has also been reviewed 
in Ref. [13], and more recently in Ref. [15] . 

The remainder of this article is structured as follows: we start by giving a more detailed discus- 
sion of proper magnetic ferroelectrics, and we summarize some of the key developments where 
first principles studies have made important contributions. We then focus in particular on re- 
search related to BiFe03, which is probably the most studied multiferroic material to date. After 
that, we give an overview over more recent advancements in the field of improper multiferroics, 
and discuss some recent work on various manganite systems: orthorhombic E-type H0M11O3 and 
half-doped Lao.sCao.sMnOs. We end with some conclusions and perspectives for future research. 

Finally, before starting our discussion of proper and improper multiferroics, we want to mention 
that even though no new calculational techniques have to be developed for the study of these 
materials, research on multiferroics typically involves a combination of a variety of advanced 
techniques, most of which have been established only during the last decade (roughly speak- 
ing). These techniques include for example beyond-LDA/GGA approaches for the treatment of 
strongly correlated transition metal oxides, mostly LSDA+C/ |161I17|. methods for the treatment 
of non-collinear magnetism [T5] and spin-orbit coupling [TS], the Berry phase approach to cal- 
culate electric polarization [20\ [21] combined with a further analysis using maximally localized 
Wannier functions [22], and many more. 

2 Proper magnetic ferroelectrics 

Most of the "early" first principles work on multiferroics was focused on proper magnetic fer- 
roelectrics, in particular on identifying mechanisms for ferroelectricity that are compatible with 
the simultaneous presence of magnetic order. 

In conventional ferroelectrics such as BaTiO"3 and PbTi03, hybridization effects between the 
filled oxygen p states and the empty transition metal d states are essential for the appearance of 
the structural instability that causes ferroelectricity [23] . Early first principles work pointed out 
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that such a mechanism is unfavorable if the transition metal d states are partially filled, which 
to some extent explains the relative scarcity of magnetic ferroelectrics \2A\ [25] . 

The ferroelectricity in multiferroic materials is therefore generally caused by a different mech- 
anism than in prototypical ferroelectric materials such as BaTi03, PbTi03, or KNbC>3, which 
all contain transition metal cations with a formal d° configuration. As in the case of these 
conventional ferroelectrics, electronic structure calculations have been crucial in identifying and 
classifying different mechanisms for ferroelectricity that are also compatible with the simultane- 
ous presence of partially filled d of / states. 

Two such mechanisms have emerged from these early studies: 

1. Ferroelectricity caused by stereo-chemically active "lone-pair" cations, e.g. Bi 3+ or Pb 2+ . 

2. "Geometric ferroelectricity", where the structural instability is driven by size effects and 
other geometrical considerations. 

It is well known in chemistry, that cations containing a highly polarizable 5s or 6s lone pair 
of valence electrons have a strong tendency to break local inversion symmetry. This can be 
understood by a mixing between ns and np electron states, which can lower the energy of the 
cation, but is only allowed if the ionic site is not an inversion center. Alternatively, in a solid 
this tendency can be understood as cross-gap hybridization between occupied oxygen p and 
unoccupied np states of the lone-pair cation, similar to the cross-gap hybridization between 
occupied oxygen p and unoccupied transition metal d states that gives rise to the ferroelectricity 
in conventional ferroelectrics [26J. In fact, the presence of the lone-pair active Pb 2+ cation is an 
important factor for the ferroelectric properties of PbTi03 (in addition to the presence of the 
d° Ti 4+ cation) [23]. The lone -pair mechanism was identified as the source of the ferroelectric 
instability in BiMn0 3 [271 EE] and BiFe0 3 [SUED]. 

In contrast to this, the ferroelectric instability in geometric ferroelectrics does not involve any 
significant re-hybridization effects. Instead, a structural instability in such systems is generated 
mainly by size effects and geometric constraints, i.e. the space-filling and ionic coordination in 
the "ideal" high-symmetry structure is not optimal, but can be improved by a small distortion 
that eventually breaks inversion symmetry. The first material that was identified as geometric 
ferroelectric is hexagonal YMn03 [31] (see Fig. [Tk) . First principles calculations showed that the 
ferroelectric structure of this material results from an interplay between a polar T-point mode and 
a non-polar Brillouin zone-boundary mode that leads to a unit cell tripling |3H [7] . Furthermore, 
calculated phonon frequencies together with group theoretical analysis suggests that YMnO"3 
is an improper ferroelectric, where the hexagonal point group of the centrosymmetric high- 
symmetry structure allows a coupling between the otherwise stable and the unstable K% 
mode [TJ. 

An example for proper geometric ferroelectricity has been found in the series of antiferromagnetic 
(AFM) fluorides BaMF4, where M can be Mn, Fe, Co, or Ni [32J. The special connectivity of 
the fluorine octahedra in these systems, which are arranged in quasi-two-dimensional sheets, 
gives rise to one unstable phonon mode that involves alternating octahedral rotations together 
with an overall shift of the interjacent Ba cations relative to the other ions (see Fig. QJ)). This 
shift creates an electric dipole moment, and since only one structural mode is involved the 
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Figure 1: Crystal structures of various magnetic ferroelectrics: a) YMn03, which has been clas- 
sified as improper geometric ferroelectric, crystallizes in a layered hexagonal structure, consist- 
ing of a two-dimensional arrangement of connected oxygen bi-pyramids surrounding the Mn 3+ 
cations that are separated by layers of Y 3+ cations, b) BaNiF4, a proper geometric ferroelec- 
tric, is found in an orthorhombic structure with buckled planes of fluorine octahedra around 
the Ni 2+ cations and additional interjacent Ba 2+ cations, c) BiFeOs, where the ferroelectricity 
is driven by the stereochemically-active Bi 3+ cation, exhibits a rhombohedrally distorted per- 
ovskite structure, where all ionic sublattices are displaced relative to each other along the polar 
(111) direction, and the oxygen octahedra are rotated around the same (111) axis, alternately 
clockwise and counter-clockwise. 



5 



corresponding ferroelectricity is classified as "proper". Due to the fact that fluorine systems 
are generally much more ionic and less covalent than oxides, geometric ferroelectricity can be 
expected to be the dominant source for ferroelectric instabilities in fluoride compounds. 

Very recently, the question of why exactly the standard p-d hybridization mechanism for ferro- 
electricity is unfavorable for systems with partially filled d shells has been revisited |33(. 134]. For 
perovskite systems, with dominantly cubic crystal field splitting between the t2 g and e g mani- 
folds, it is not fully clear why for example a d 3 configuration with partially filled t2 g states, but 
empty e g orbitals, cannot give rise to a favorable cross-gap hybridization between filled oxygen 
p and empty transition metal e g states. It was suggested that the Hund's coupling between t2 g 
and e g states will disfavor such hybridization [11] . This was supported by LDA+t/ calculations 
for CaMn03, where the Hund's coupling was effectively "turned off", which indeed resulted in a 
tendency for off-centering of the Mn 4+ cation. In addition, recent first principles calculations for 
CaMn03, SrMn03, and BaMn03 in the perovskite structure show that these systems can de- 
velop a ferroelectric instability, but that this ferroelectric instability competes with a non-polar 
"antiferrodistortive" instability, and that the relative strength of these two instabilities depends 
strongly on the unit cell volume |33[ I34j. For larger volumes (i.e. BaMnOs) the ferroelectric 
instability becomes dominant. Thus, even though BaMn03 is not stable in the cubic (or in 
the orthorhombically distorted) perovskite structure (it crystallizes in a hexagonal structure), 
this opens up the possibility to stabilize the corresponding ferroelectric phase by using epitaxial 
constraints, i.e. using thin film growth techniques. 

Apart from these investigations into possible mechanisms for ferroelectricity that are compatible 
with the simultaneous presence of magnetic order, first principles calculations have also been used 
to rationalize experimental observations, investigate possible mechanisms for coupling between 
the electric polarization and the magnetic order, and to design new multiferroic and magneto- 
electric materials. In the following we will highlight some of these calculations, in particular the 
work related to one of the most prominent multiferroic materials: bismuth ferrite. 

2.1 First principles calculations for BiFeC>3 and related work 

BiFeC>3 (BFO) is one of the most studied (probably the most studied) multiferroic material. 
BFO is known to be multiferroic (or more precisely: AFM and ferroelectric) already since the 
early 1960s [35]. However, for a long time it was not considered as a very promising material for 
applications, since the electric polarization was believed to be rather small [36] and the AFM 
order does not lead to a net magnetization |37[ [38] . 

This has changed drastically, following a publication in Science in 2003 (Ref. [39]), which to great 
extent has triggered the intensive experimental and theoretical/computational research on BFO 
during the last 5-6 years. In this study, a large spontaneous electric polarization in combination 
with a substantial magnetization was observed above room temperature in thin films of BFO 
grown epitaxially on SrTi03 substrates. The presence of both magnetism and ferroelectricity 
above room temperature, together with potential coupling between the two order parameters, 
makes BFO the prime candidate for device applications based on multiferroic materials. 

Whereas the large electric polarization was later confirmed independently, and explained by first 
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principles calculations, the origin of the strong magnetization reported in [39] is still unclear 
and, to the best of our knowledge, it has never been reproduced in an independent study. It 
is generally assumed that the magnetization reported in Ref. [39] is related to extrinsic effects 
such as defects or small amounts of impurity phases. 

The large electric polarization, which appeared to be at odds with bulk single crystal measure- 
ments from 1970 [36], was originally assumed to be due to epitaxial strain, which results from 
the lattice constant mismatch between BFO and the substrate material SrTi03. It is known 
that epitaxial strain can have drastic effects on the properties of thin film ferroelectrics. For 
example, it can lead to a substantial enhancement of electric polarization and can even induce 
ferroelectricity at room temperature in otherwise non- ferroelectric SrTiOs |40|, I41j. 

In the following we illustrate how first principles calculations have been instrumental in clarifying 
the origin of both polarization and magnetization in thin film BFO, by showing that the large 
electric polarization found in the thin films is in fact intrinsic to unstrained bulk BFO and that, 
in contrast to many other ferroelectrics, epitaxial strain has only a minor effect in this material. 

2.1.1 Electric polarization of bulk BFO and the effect of epitaxial strain 

According to the so-called "Modern theory of polarization", the electric polarization of a bulk 
periodic system is defined via the Berry phase of the corresponding wavefunctions [2U1 I2"T] . 
Since this geometrical phase is only well defined modulo 2tt, the polarization is only well-defined 
modulo so-called "polarization quanta", given by AP = jfdi, where e is the electronic charge, 
Si a primitive lattice vector (i = 1, 2, 3), f2 the unit cell volume, and / is a spin degeneracy factor 
(/ = 2 for a non-spinpolarized system, / = 1 for a spin-polarized system). If the expression 
for the polarization is recast as a sum over "Wannier centers" [20], a translation of one of the 
occupied Wannier states from one unit cell to the next corresponds to a change in polarization by 
exactly one "quantum" . The multivaluedness thus reflects the arbitrary choice of basis vectors 
when describing an infinite periodic structure. 

In spite of this multivaluedness of the bare polarization for a specific atomic configuration, 
differences in polarization are well defined quantities, provided the corresponding configurations 
can be transformed into each other in a continuous way and the system remains insulating along 
the entire "transformation path" |21j . 

In particular, the spontaneous polarization of a ferroelectric material is defined as half the differ- 
ence in polarization between two oppositely polarized states, or equivalently, as the difference in 
polarization between the ferroelectric structure and a suitable centrosymmetric reference config- 
uration. In order to calculate the spontaneous polarization one therefore has to perform a series 
of calculations for different configuration between the ferroelectric state and the centrosymmet- 
ric reference structure. If the change in polarization between two such configurations is much 
smaller than the polarization quantum, then the corresponding difference can be clearly iden- 
tified and the full change in polarization along the transformation path, i.e. the spontaneous 
polarization, can be determined. 

The application of this procedure to calculate the spontaneous polarization of BFO is compli- 
cated by the following two features: i) the polarization quantum for a spin-polarized system is 
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Figure 2: Evolution of the polarization P along the transformation path from a negatively 
polarized state ( — 100 % distortion), through a centrosymmetric reference configuration (0 % 
distortion), to a positively polarized state (+100 % distortion). Red circles correspond to the 
LSDA+C/ calculation with f7 e g = 2 eV, green diamonds indicate the LSDA result for the fully 
polarized states. Different values of P for fixed amount of distortion are separated by the 
polarization quantum APq 111 ^ = 186 /xC/cm 2 . The spontaneous polarization P s is given by the 
difference in polarization between the fully distorted and the undistorted configuration for an 
arbitrary branch of the bare polarization. Note: the systematic sketches at the bottom do not 
correspond to the actual crystal structure of BFO. 



only half that for a similar nonmagnetic system, and ii ) due to the underestimation of the local 
spin splitting for Mott-Hubbard insulators within the standard local spin-density approximation 
(LSDA), BFO becomes metallic for the less distorted reference configurations within LSDA. 

These problems have been overcome in Ref. [29] by using the LSDA+J7 method [T7] to 
calculate the electronic structure of BFO in various configurations along the transformation 
path from the fully distorted R3c structure to the centrosymmetric cubic perovskite (Pm3m) 
structure. Within the LSDA+J7 method the local d-d exchange splitting is enhanced by the 
Hubbard U and BFO stays insulating even in the undistorted cubic perovskite structure (for U 
values U eS = U-J = 2-4 eV [29]). 

Figure [2] shows the evolution of the electric polarization with varying degree of distortion between 
two oppositely polarized states calculated for U c r = 2 eV. The LSDA results are included for 
±100 % distortion. The fact that the corresponding symbols (green diamonds) can barely be 
recognized behind the red circles that indicate the LSDA+C7 results shows that the value of the 
bare polarization is rather insensitive to the exact value of U c q. It can be seen that different 
values of P corresponding to the same amount of distortion are separated by the polarization 



8 




8 [%] 

Figure 3: Dependence of the spontaneous polarization P s on epitaxial strain e for BFO in two 
different structural modifications and some other (non-magnetic) ferroelectrics. Symbols corre- 
spond to results from first principles calculations for strained unit cells (data for BaTiOa/PbTiOa 
is taken from [13] /[S]), lines are obtained from the calculated bulk linear response functions 
(see |45|). Note that the epitaxial constraint for all systems is applied in the plane perpendicular 
to the polarization, i.e. (001) for BaTi0 3 , PbTi0 3 , and P4mm-BiFe0 3 , and (111) for LiNb0 3 
and P3c-BiFe0 3 . 

quantum along (111), APq 111 ^ = f^(a\ + S2 + S3), where 01,2,3 are the primitive lattice vectors of 
the rhombohedral R3c structure. As indicated, the spontaneous polarization P s can be obtained 
as the difference between the fully distorted and the undistorted configuration for an arbitrary 
"branch" of the bare polarization. 

From these calculation a spontaneous polarization of bulk BFO of ~ 95 //C/cm 2 has been 
obtained. This is an order of magnitude larger than what was previously believed to be the 
case, based on the measurements in Ref. [36], and even exceeds the polarization of typical 
prototype ferroelectrics such as BaTi0 3 , PbTi0 3 , or PbZro.5Tio.5O5 (PZT). Variation of U e g 
within reasonable limits changes the calculated value for the electric polarization by only ~ 
±5 /iC/cm 2 , i.e. the large value of the polarization is rather independent from the precise value 
of the Hubbard parameter. This is consistent with the assumption that the transition metal d 
states do not play an active role for the ferroelectric instability in BFO. The calculated large 
spontaneous polarization for bulk BFO is also consistent with the large ionic displacements in 
the experimentally observed R3c structure of BFO (see Fig. [Tb), compared to an appropriate 
centrosymmetric reference configuration. Recently, the large polarization of ~ 100 fiC/cm 2 
along (111) for bulk BFO has also been confirmed experimentally by new measurements on 
high-quality single crystals |42j . 
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Effects of epitaxial strain can be assessed from first principles by performing bulk calculations 
for a strained unit cell, where the lattice constant within a certain lattice plane (corresponding 
to the orientation of the substrate surface) is constrained, whereas the lattice constant in the 
perpendicular direction as well as all internal structural parameters are allowed to relax. Such 
calculations have been performed for BFO corresponding to a (111) orientation of the substrate 
[46j . In this case the R3c symmetry of the bulk structure is conserved and the epitaxial con- 
straint is applied in the lattice plane perpendicular to the polarization direction. It was found 
that the sensitivity of the electric polarization to strain is surprisingly weak in BFO, much 
weaker than in other well-known ferroelectrics [46] (see Fig. [3]) . A systematic comparison of the 
strain dependence in various ferroelectrics, including BFO in both the R3c and a hypothetical 
tetragonal phase with PAmm symmetry, has been performed in Ref. [35] (see Fig. [3]). It was 
shown that the effect of epitaxial strain for all investigated systems can be understood in terms 
of the usual bulk linear response functions and that both strong and weak strain dependence 
can occur. 

Systematic calculations corresponding to a (001) orientation of the substrate, the one that is 
most often used experimentally, have not been performed so far. Since the epitaxial constraint 
in this case breaks the rhombohedral symmetry of the bulk structure, the corresponding strained 
unit cell has a lower symmetry with more free parameters than in the (lll)-strained case. Nev- 
ertheless, the effect of such a monoclinic strain on the ferroelectric polarization in BFO has been 
investigated by performing calculations for a set of lattice parameters derived from represen- 
tative experimental data. Due to the lower symmetry, the polarization in this case is slightly 
rotated away from the (111) direction, but the overall magnitude remains nearly unchanged 
compared to the unstrained case. From this it was concluded that the polarization in BFO 
is generally rather insensitive to epitaxial strain, and that the large polarization measured in 
thin films is basically the same as in the corresponding bulk system. Indeed, the polarization 
of ~ 60 /iC/cm 2 reported in Ref. [39j for a (001) oriented thin film agrees well with the corre- 
sponding projection of the calculated bulk value (which is oriented along the (111) direction), 
and polarization measurements for BFO films with different substrate orientations ((001), (101), 
and (111)) can all be understood by assuming that the polarization vector in all cases points 
essentially along (111) and has approximately the same length [57]. More recently, systematic 
experimental investigations of the strain effect in epitaxial BFO films have been undertaken by 
comparing results of BFO films with different thicknesses, which have confirmed the predicted 
weak strain dependence of the polarization in BFO [38] 

Finally, it should be noted that Ref. [39] also contains results of first principles calculations for the 
electric polarization of two structural variants of BFO: the rhombohedral bulk structure with R3c 
space group, and a hypothetical tetragonal structure with PAmm symmetry, based on the lattice 
parameters found in the thin film samples. At that time it was assumed that such a tetragonal 
phase is stabilized in epitaxial thin films and that the difference in polarization observed in 
thin films compared to bulk BFO was due to a large difference in polarization between the two 
different structural modifications. However, the DFT results presented in Ref. [39J were not 
conclusive, since only the bare polarization for the two different structures was reported and not 
the spontaneous polarization that is measured in the corresponding "current-voltage" switching 
experiments. 
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In fact, it is indeed possible that a different phase is stabilized in thin films, which can then 
lead to more significant changes of ferroelectric and magnetic properties compared to bulk BFO. 
However, it is important to distinguish between the simple case of a somewhat distorted version 
of the rhombohedral bulk structure and a truly different phase, which would for example be 
characterized by a different oxygen octahedra tilt pattern or a different number of formula units 
contained in the crystallographic unit cell. 

Calculations presented in Ref. [35] (see also [M]) show that if BFO is constrained to tetragonal 
PAmm symmetry (with no octahedral tilts and only one formula unit per unit cell) it develops a 
" super-tetragonality" with c/a ratio of 1.27 and a giant electric polarization of P s ~ 150 ^C/cm 2 . 
A polarization of this magnitude has indeed been found in some highly strained films with c/a 
ratios between 1.2-1.3 |50|. 149] . whereas many other experimental reports of "tetragonal" BFO 
films with smaller c/a ratio also exist. These reports should be regarded with some caution, since 
the structural characterization of thin films is usually restricted to the measurement of lattice 
constants and of angles between certain crystallographic directions. A full characterization 
of ionic distortions (including octahedral tilt patterns etc.) is generally not possible for thin 
films, and first principles calculations can therefore play an important role in clarifying open 
questions about the exact thin film structure of BFO. In principle, if one tries to epitaxially 
match the rhombohedral bulk structure of BFO on a square lattice substrate plane, one can 
expect to obtain a mono clinically distorted version of the BFO bulk structure. However, since 
the rhombohedral angle in bulk BFO is very close to 60°, the value that corresponds to an 
underlying cubic lattice, the monoclinic distortion can be rather small, and the thin films might 
appear essentially tetragonal. 

2.1.2 Weak ferromagnetism in thin film BFO and coupling between the various 
order parameters 

In addition to these structural studies, DFT calculations have also been used to investigate the 
magnetic properties of BFO, in particular the possible origin for the significant magnetization 
reported in Ref. [Hj5]. Bulk BFO is known to exhibit "G-type" AFM ordering [37], i.e. the mag- 
netic moment of each Fe cation is antiparallel to that of its nearest neighbors. Superimposed to 
this G-type magnetic order a long-period cycloidal modulation is observed, where the AFM or- 
der parameter L = M1 — M2, defined as the difference between the two sublattice magnetizations 
Mi 5 2, rotates within the (110) plane with a wavelength of ~ 620 A [38] . 

Calculations for bulk BFO show a very strong and dominant AFM nearest neighbor interaction 
[51] , in agreement with the observed G-type magnetic order and the rather high Neel temperature 
of oj 600 K. In addition, the magnetocrystalline anisotropy has been calculated, and a preferred 
orientation of the Fe magnetic moments perpendicular to the polar [111] direction has been 
found [52]. Within the (111) plane a 12-fold degeneracy remains, leading to an effective "easy- 
plane" geometry for the magnetic moments. For an orientation of the AFM order parameter 
L within this (111) plane, weak ferromagnetism is symmetry-allowed, i.e. a small canting of 
the two AFM sublattice magnetizations can occur, which results in a net magnetization [53] . 
Indeed, if spin-orbit coupling is included in the calculation (while the cycloidal modulation is 
neglected), a small canting of the magnetic moments is obtained [52]. The magnitude of the 
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Figure 4: Dependence of the weak magnetization in BFO on the LSDA+C7 parameters U c s = 
U — J and J. The dash-dotted line represents the reported value of 0.05 /iB/Fe. The sketches on 
the left side illustrate how the canting of the two AFM sublattice magnetizations, represented 
by the magnetic moments Mpei and Mp e 2 of the two Fe cations in the primitive unit cell, gives 
rise to the net magnetization M. 

resulting magnetization depends on the choice of the Hubbard U and the Hund's rule parameter 
J, but for reasonable values of U e s = U — J the magnetization is around 0.05 /is/Fe cation (see 
Fig. [3D . This value of the magnetization agrees quite well with various thin film measurements 
|54[ [55"1 [56] , but is significantly smaller than what was originally reported in Ref . [39] . It has 
to be pointed out that no magnetization is observed in bulk BFO, where the presence of the 
cycloidal modulation effectively cancels any net magnetic moment. If the cycloidal modulation 
is suppressed, either by applying a strong magnetic field [57] or by chemical substitution [SS] 
a small magnetization appears, with comparable magnitude to the computational result. It is 
generally assumed that the cycloidal rotation of the AFM order parameter is also suppressed 
in thin films, likely due to enhanced anisotropy, and that the small magnetization observed in 
the thin films is due to weak ferromagnetism. This is supported by a neutron diffraction study 
on BFO films, which could not find the satellite peaks associated with the cycloidal modulation 

m- 

Furthermore, first principles studies addressing the effect of epitaxial strain and the presence 
of oxygen vacancies did not find a significant increase in magnetization [35], and it is therefore 
likely that the large magnetization reported in [39] is due to other defects or small amounts of 
impurity phases. 

The appearance of weak ferromagnetism in thin films of BFO leads to the question of whether 
this small magnetization is coupled to the electric polarization, i.e. whether it can be manipu- 
lated by applying external electric fields. Indeed, the absence of an inversion center located at 
the midpoint between two interacting magnetic moments is crucial to produce a non-vanishing 
Dzyaloshinskii-Moriya (DM) interaction, which has been identified as the microscopic mecha- 
nism responsible for the magnetic moment canting in weak ferromagnets [59]. Thus, inversion 
symmetry breaking can cause both weak ferromagnetism and ferroelectricity, suggesting possible 
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cross-correlations between these two properties. First principles calculations have been used to 
explore this possibility for magnetization-polarization coupling in BFO [52] and in BaNiF4 |60j . 
It was found that in BFO the DM interaction is caused by a non-polar antiferrodistortive mode, 
not by the polar distortion, and therefore the weak ferromagnetism in BFO is not controlled by 
the spontaneous polarization and cannot be switched using an electric field [52]. In contrast, in 
BaNiF4, it is indeed the polar distortion that creates a DM interaction, but the symmetry is 
such that no net magnetization results. Instead, a secondary (weak) AFM order parameter is 
induced in addition to the distinctly different primary AFM order [60]. Only recently, a mate- 
rial has been suggested, based on a combination of first principles calculations and symmetry 
considerations, that fulfills all requirements for "ferroelectrically-induced weak ferromagnetism" 
[61j . The corresponding material, R3c structured FeTiOs, is closely related to BFO in that 
it has the same overall structural symmetry, but with the magnetic Fe cations located on the 
perovskite A site instead of the perovskite B site as in BFO. It is this difference in the local 
site symmetry of the magnetic cation, that is crucial for the coupling between the spontaneous 
polarization and the weak magnetization |6H 162] . Experimental work is currently underway to 
validate this theoretical prediction. 

2.1.3 Designing new multiferroics and new functionalities 

The prediction of FeTi03 as a possible candidate for electric field switchable weak ferromag- 
netism, is one example for attempts to design new materials with novel or more favorable 
magneto-electric properties based on first principles electronic structure calculations. 

Another example is the design of a material that allows for magneto-electric phase control [63]. 
Calculations for the rare-earth magnet EuTiOs showed that this material exhibits a soft infrared- 
active, i.e. polar, phonon mode that becomes unstable if the material is epitaxially strained. In 
addition, due to strong spin-phonon coupling in this material, the instability is more pronounced 
for ferromagnetic ordering of the Eu spins than for the case of an AFM arrangement. Since the 
ground state magnetic structure for the lower strain region is AFM, it was suggested that a 
phase transition from a non-polar AFM phase into a ferroelectric-ferromagnetic phase can be 
induced by applying a strong magnetic field, if the material can be prepared in thin films with 
a compressive epitaxial strain of around 1 % [63] . 

In addition, attempts have been made to design materials that combine strong ferroelectric 
polarization with a large magnetization above room temperature. If such a material would 
also exhibit pronounced coupling effects between polarization and magnetization, which ideally 
would allow to switch the polarization via a magnetic field or vice versa, then this would probably 
create a similar excitement as finding a room temperature superconductor. Unfortunately, at 
the moment no multiferroic that exhibits all these properties is known (similarly, no room 
temperature superconductor is known at present). 

A suggestion for a material combining large polarization and large magnetization has been 
made in Ref. [63]. First principles calculations predict, that if half of the Fe 3+ cations in BFO 
are replaced by Cr 3+ cations in a checkerboard-like ordered arrangement, then the resulting 
material Bi2FeCr06 is stable in a rhombohedral structure similar to BFO with a spontaneous 
ferroelectric polarization of around 80 jitC/cm 2 and a magnetization of 2 fi^ per formula unit. 
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The magnetization in this case results from a ferri-magnetic arrangement, where the magnetic 
moments of the Cr cations are antiparallel to those of the Fe cations. A subsequent study of 
the strength of the magnetic coupling in the series of compounds BiFe03-Bi2FeCr06-BiCr03 
has found that the Neel-temperature in Bi2FeCrOg is unlikely to be above room temperature 
[51j . but nevertheless several attempts have been made to synthesize the corresponding material 
[651 166} 167]. The synthetic challenge here, is to achieve the required checkerboard-type ordering 
of Fe and Cr cations on the B sites of the underlying perovskite structure, which might be 
possible by utilizing layer-by-layer growth on a (lll)-oriented substrate. 

2.2 Perspectives for future studies of proper multiferroics 

The examples discussed so far show that first principles calculations have proven not only to 
be useful for rationalizing experimental observations and identifying different mechanisms for 
ferroelectricity that can be found in multiferroic materials, but also to facilitate quantitative pre- 
dictions of new materials and novel effects in proper magnetic ferroelectrics. Future applications 
of ab initio methods in the design of new materials and in calculating the expected properties 
of these materials are therefore expected to continue to have a significant impact on the overall 
progress of this field. 

In particular, a material with large magnetization and large polarization above room temperature 
is still elusive. From the current point of view there is no fundamental reason why such a material 
should not exist, and creative ideas on how to circumvent the limitations and restrictions of 
materials chemistry that have been encountered so far are still highly desirable. 

Another area where DFT will undoubtedly have (and already has) a substantial impact, is the 
study of artificial heterostructures consisting of a combination of magnetic and ferroelectric 
materials [13] , Examples of computational work in that direction that have already appeared 
include the study of artificial tri-layered superlattices of different magnetic and nonmagnetic 
oxides [US] and the investigation of polarization effects at the interface between a ferromagnetic 
metal and a ferroelectric insulator [69]. 

3 Improper Multiferroics 

In the beginning of this section, we will focus on the origin of ferroelectricity in the so-called 
"Improper multiferroics" (IMF), outlining a few differences with respect to the more conventional 
"proper" multiferroics discussed so far. 

As pointed out in the previous sections, in displacive ferroelectric materials (such as prototypi- 
cal perovskite-like BaTi03 or multiferroic BiFeOa), due to strong covalency effects, the relative 
displacement of the anionic sublattice with respect to the cationic sublattice gives rise to a spon- 
taneous and switchable polarization, which is the (primary) order parameter in the ferroelectric 
transition. On the other hand, in IMF, the primary order parameter of the phase transition 
is related to electronic (i.e. spin, charge, or orbital) degrees of freedom [12] . The important 
thing is that the resulting electronic order lacks inversion symmetry (IS), therefore opening the 
way to ferroelectricity. Therefore, polarization occurs as a by-product of the electronic phase 
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transition and can be described as a "secondary" order parameter. As a consequence, i) even 
the state with ions pinned in centrosymmetric positions can show a finite (purely electronic) 
polarization; ii) the ions can "react" to the non-centrosymmetric charge-redistribution by dis- 
placing, so as to give a (more traditional) ionic contribution to the total polarization. In order 
to push ahead with the comparison between proper and improper multiferroics, one can say that 
ferroelectricity in IMF is driven by "correlation" effects (as related to spin or charge arrange- 
ments), at variance with the previously mentioned case of standard ferroelectrics where it is 
mostly driven by covalency. In IMF where polarization is magnetically-induced, it is reasonable 
to expect a strong coupling between magnetic and ferroelectric properties, since the two dipolar 
and magnetic orderings share the same origin and occur at the same temperature. 

In Fig. we schematically classify IMF on the basis of the different mechanisms to induce 
ferroelectricity that have been proposed so far. We would like to point out that what we present 
in the following is a non-exhaustive list of the IMF materials and related mechanisms. In fact, 
IMF represent a quickly evolving field: new materials and/or novel mechanisms are proposed 
on a monthly or even weekly basis. With no doubt, we therefore expect in the near future this 
classification to become richer in compounds and to expand as far as mechanisms are concerned. 

In Fig.OlMF are divided in two main classes: those where ferroelectricity is driven by spin-order 
(i.e. where the "magnetic" arrangement breaks IS) and those where it is driven by charge-order 
(i.e. where the charge-disproportionation leads to a non-centrosymmetric arrangement). In 
turn, the magnetically-induced ferroelectricity can occur in two different ways: i) the first and 
most studied case where a non-collinear spin-spiral occurs and the IS-breaking arises due to 
a spin-orbit related mechanism in the DM-like antisymmetric exchange term [70 1 P7T| I72j : ii) 
the case of (mostly collinear) AFM spins where the IS-breaking occurs in the Heisenb erg-like 
symmetric exchange-term [73 1 [73] . 

Along with the classification of IMF, we show in Fig. [5] a few links to IMF materials for which 
ab-initio studies have been reported in the literature. 

Chronologically, the recent interests towards IMF were boosted by the discovery of ferroelectric- 
ity in TbMn03 and of the control of the polarization direction achieved via an applied magnetic 
field [87]. However, the ab-initio simulations for TbMn03 came much later |76| I75|. due to the 
complexity in the related simulations: advanced capabilities (such as non-collinear magnetism 
and spin-orbit coupling) are needed to reproduce the observed tiny effects, which implicitly re- 
quires a high precision in terms of numerical parameters in the calculations. In the ab-initio 
field, the first IMF to be studied were collinear antiferromagnets, such as TbM^Os [ST] and 
HoMn03 [74] . Since the latter will be described in detail in Sec. 13.1.11 we will now briefly 
discuss the first one. The class of manganites often labeled as "1-2-5" from the stoichiometry 
of rare-earth, transition metal, and oxygen, respectively, is an actively studied set of IMF. De- 
spite some non-collinearity and non-commensurability effects, most of the mechanisms behind 
multiferroicity can be described through simulations with non-centrosymmetric collinear spin 
arrangement using a relatively small supercell. The suggested polarization was of the order of 
1 f-iC/cm 2 and the polarization was reversed by changing the spin-orientation in the unit cell, 
providing evidence for the magnetic origin of ferroelectricity in TbMi^Os. Within the same 
class of materials, HoM^Os was studied in Ref. [82]: the main and new result of that work was 
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Figure 5: Schematic classification of IMF, in terms of different mechanisms (left side) and 
compounds (right side). The (non comprehensive) list includes a few materials which were 
studied by first-principles (see related references: a) Ref. [75l [76], b) Ref. [77], c) Ref. [741 [78] . 
d) Ref. [7J[80], e) Ref. [ED [82], f) Ref. [83], g) Ref. [M], h) Ref. [85], i) Ref. [86]). 



that the ionic and electronic contributions were strongly dependent on the value of the Hubbard 
U parameter used in a LSDA+C7 approach, pointing to the important role of correlation effects 
in 1-2-5 manganites. 

Within the spin-spiral class of IMF, Li-Copper-based oxides were the first compounds to be 
studied from first-principles [77]: upon switching-on spin-orbit coupling, the calculated polar- 
ization was rather small (of the order of tens or hundreds of //C/cm 2 , depending on whether 
ionic relaxations were included or not in the simulations). Shortly later, the prototypical case of 
TbMnOs was published in two important papers (one following the other in Phys. Rev. Lett.), 
Refs. |76t 175]. It was shown that the purely electronic contribution (i.e. evaluated by switching 
on spin-orbit but keeping the ions frozen into their paramagnetic centrosymmetric configura- 
tion) was much smaller than the ionic contribution (i.e. evaluated by relaxing the ions). In the 
TbMnC>3 case, the order of magnitude of the ab-initio polarization was found to be in excellent 
agreement with experiments [87]. Remarkably, at the time of publication, the sign of polariza- 
tion obtained within DFT was opposite with respect to experiments; indeed, it later turned out 
|88| that the discrepancy was due to a misunderstanding in the conventions of the experimental 
settings and an excellent agreement between theory and experiments could be finally obtained. 

Within the field of charge-order-induced ferroelectricity, a prototype has emerged: the triangular 
mixed- valence iron-oxide, LuFe204 |89| . There, the frustrated charge-ordering is such as to lack 
centrosymmetry: in each FeO bilayer, there is an alternation of iron atoms, with Fe 2+ :Fe 3+ 
ratios of 2:1 and 1:2, therefore giving rise to a polarization within each bilayer. The polarization 
estimated from first-principles is very large (of the order of 10 fiC/cm 2 in the bilayer). However, 
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some controversy exists for that material, since it is questioned whether the stacking of the 
bilayers is such as to produce net ferroelectricity [83] or a global antiferroelectricity with no net 
polarization [90]. More work (both from theory and from experiments) will be needed in that 
respect. 

Recently, another collinear compound has been studied, CasCoMnOg [7SJ[80]. The main ab-initio 
findings were: i) a large Co orbital moment, which renders the system similar to an Ising-like 
chain, with alternating trigonal prismatic Co 2+ and octahedral Mn 4+ sites in the spin chain; ii) 
a large calculated polarization (about 1.7 /iC/cm 2 ), caused by a significant exchange-striction 
combined with a peculiar TTJ-J- spin configuration. 

Given this general background, in the following sections we will present some examples of ab- 
initio calculations for IMF. In closer detail, we will discuss rare-earth manganites (c/r. Sec. 13. 1.1]) 
[7H [78] and hole-doped manganites (cfr. Sec. 13.1 . 2j) [85] as examples of AFM materials where 
the spin-arrangements break inversion symmetry, with polarization being due to Heisenb erg-like 
mechanisms. We will conclude the section by discussing some perspectives and open issues in 
the field. 

In what follows, we will mainly show the results of DFT simulations performed using the Vienna 
Ab-initio Simulation Package (VASP) [ST] and the generalized gradient approximation [92] to the 
exchange-correlation potential. For the construction of the Wannier functions, we used the Full- 
potential Linearized Augmented Plane- Wave (FLAPW) [S3] code in the FLEUR implementation 
[94]. For a better treatment of correlation effects, the so-called LSDA+£7 approach [17] (with 
U = 4 eV and J = 0.9 eV) was used in the case of hole-doped manganites. For further technical 
details, as far as computational or structural parameters are concerned, we refer to our original 
publications [TJEEIES]. 

3.1 Highlights on Improper Multiferroics 
3.1.1 E-type rare-earth ort ho- manganites 

Let us start the discussion of ferroelectricity in orthorhombic manganites, -RMnOs, by plotting 
the AFM spin-arrangement characteristic of the E-type HoMn03 . In Fig. [B^, we sketch the 
ions in the MnO"2 plane and highlight the zig-zag spin-chains, typical features of the E-type 
antiferromagnetism: zig-zag ferromagnetic (FM) spin-up-chains (green atoms in Fig. [6^) are 
antiferromagnetically coupled to neighboring spin-down-chains (pink atoms in Fig. [6^,) . The out- 
of-plane coupling is also AFM. We note that the antiferromagnetically-coupled zig-zag chains 
lead to a doubling of the conventional GdFe03-like unit cell (20 atoms, Pnma space group) 
along the a-axis. Indeed, the E-type was experimentally observed to be the magnetic ground 
state in distorted manganites with small ionic radius for the rare-earth ion (i.e. R = Ho, . . . , Lu) 
|96t I97j. It was shown |12l 178] that the stabilization of an TTII spin-chain (as the one present in 
the E-type along the diagonal directions in the a-c plane, cfr. Fig. [Bfci), is driven by i) a relatively 
small nearest-neighbor exchange coupling constant; ii) a large AFM next-nearest-neighbor; Hi) 
a quite large magnetic anisotropy so that the spins can be considered as Ising-like. 

Why should the E-type magnetic configuration lead to a ferroelectric polarization? This can be 
rationalized in different (though somewhat inter-connected) ways, depending on the orbitals or 
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Figure 6: a) Ionic arrangement of AFM-E H0M11O3 in the Mn02 plane. Green (pink) rhombi 
denote in-plane projections of Mn0"6 octahedra around the up-spin (down-spin) Mn ion. Spin 
directions indicated by black arrows, b) Schematic orbital-ordering for Mn e g states. Circular 
arrows show hopping paths, as induced by the AFM-E spin configuration; green and pink arrows 
denote asymmetric hoppings for up-spin and down-spin electrons, respectively, c) Schematic 
local dipoles (denoted by blue arrows) drawn from O ap (bonded to Mn with antiparallel spins) 
to O p (bonded to Mn with parallel spins). In b) and c), the direction of polarization is also 
shown. 
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atoms one focuses on. 

Let's start with Mn e g states. Being Mn in a d 4 electronic configuration, the strong Jahn- Teller 
effect leads to two large and two small in-plane Mn-0 bond lengths, along with a staggered 
(3x 2 — r 2 )/(3y 2 — r 2 ) orbital-ordering, typical for the class of rare-earth manganites. Within a 
double-exchange-like picture, this peculiar orbital-ordering (00) leads to a favored hopping of 
the electron on the two (out of four nearest neighbors) Mn-sites towards which the orbital is 
pointing. What is peculiar of the E-type (and different from the conventional A-type in early- 
rare-earth manganites) is that, out of these two Mn sites, hopping will preferentially occur on 
the Mn with the spin parallel to the starting site, and not on the other which shows an opposite 
spin. This "asymmetric" hopping creates a "one-way path" for the electron, schematically shown 
by the circular arrows in Fig. [6b. At this point, it is clear that the short c axis is a "preferential" 
direction for the electron, with a well-defined sign for the electron hopping. This mechanism 
therefore breaks inversion symmetry and opens the way to a ferroelectric polarization P c . 

Another way to explain the direction of polarization is to look at oxygen sites. Again due to 
the peculiar E-type spin-configuration, there will be two kinds of O sites: those bonded to Mn 
with parallel spins (labeled as O p ) and those bonded to Mn with antiparallel spins (labeled as 
O ap ). Due to this inequivalency, their electronic structure will be different (even if the ions are 
frozen into a centrosymmetric "paramagnetic" configuration). This leads to a sort of oxygen 
"charge-density wave" which can be thought of in terms of a set of ordered dipoles resulting in 
a net ferroelectric component, again only along the short c-axis (c/r. Fig. [6b). 

We would now like to make one comment related to ferroelectric switching in IMF. As is well 
known, in conventional displacive perovskite-like ferroelectrics, the switched state (i.e. the 
one with —P) is achieved by displacing the ions (with respect to a reference centrosymmetric 
structure) in the opposite way compared to the +P state. However, when asking how to switch 
P in the case of magnetically-driven ferroelectrics, one might guess that some changes in the 
spin-arrangement (rather than in the ionic arrangement) should be involved. Indeed, from both 
Fig. [Bb and c, it is clear that P is switched by changing the direction of half of the spins in the 
unit cell. For example, if we revert the sign of the two spins in the central part of the unit cell 
(labelled as MnJ and Mng in Fig. [6b,), then the circular arrows in Fig. [6b will run in the opposite 
— c direction; similarly, the O-related dipoles of Fig. [6b will also change their sign. 

So far, we have taken into account purely "electronic" mechanisms, occurring when considering 
the ions frozen into their centrosymmetric configuration. However, it is reasonable to expect 
some ionic relaxations consistent with the imposed E-type spin arrangement. For example, 
according to a Heisenberg-like magnetostrictive effect, one expects that O p will try to move so 
as to gain a "double-exchange" -like energy by maximizing the Mn-O-Mn angle (recall that the 
energy lowering due to double-exchange is optimal in the ideal 180-degree case), compared to O ap 
where double-exchange is not relevant. These ionic relaxations break the atomic centrosymmetry 
and lead to an "ionic" contribution to the total ferroelectric polarization, to be added to the 
purely electronic one. 

On the basis of this introductory background, the interpretation of DFT results for HoMn03 is 
quite straightforward. It is however very important to remind that, at variance with model- 
Hamiltonian studies allowing the qualitative prediction of a selected phenomenon, first-principles 
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calculations can provide a quantitative estimate as well. Moreover, multiferroics are very com- 
plex materials where several competing mechanisms can occur. As such, identifying the strong 
and prevailing effects can be difficult within a Hamiltonian-modelling approach; on the other 
hand, all the different mechanisms are taken into account on the same footing within DFT. 

We report in Table [T] the relevant properties calculated within DFT, such as: i) the Mn-O-Mn 
angles between parallel (a p ) and antiparallel (a ap ) Mn spins, obtained after ionic relaxations 
in the presence of the E-type spin arrangement; ii) the values of the polarization calculated 
in several different ways: a purely electronic contribution (P^), estimated via the Berry- 
phase approach, when the ions are clamped in a centrosymmetric Pnma configuration; the 
polarization calculated from the so-called "Point Charge Model" (P^ M ), with the ions relaxed 
in the ferroelectric configuration, using "nominal" ionic values for the charges (i.e. 3+ on 
Mn and Ho and 2— on the O); the total (ionic + electronic) polarization in the relaxed ionic 
arrangement, calculated according to the Berry-phase approach (P^f); Hi) the Born effective 
charges, i.e. the (3,3) components of the Z* tensor for some relevant atoms: Z*(Mn), Z*(O p ) 
and Z*(O ap ). We recall that the elements are estimated by displacing the selected ion along 
the c direction by a small amount (typically about 0.01 A or less) and evaluating the change in 
the Berry-phase polarization along the same c axis. 

When focusing on the Mn-O-Mn angles, we indeed note that the angle between Mn with parallel 
spins is much larger than that where spins are antiparallel, reflecting the efficiency of relaxations 
driven by double-exchange mechanisms. As for polarization, several remarks are in order: i) 
one might naively expect a magnetically-induced mechanism to be "weak". However, this is 
contradicted by the purely electronic polarization, which is noticeably large. Moreover, this 
is one order of magnitude bigger than what was estimated in the case of spin-spirals (<~ 0.1 
/iC/cm 2 ): this reflects the efficiency of the Heisenberg vs. DM term in breaking inversion 
symmetry, ii) A similar consideration holds for the total polarization. Exchange-strictive effects 
due to the symmetric Heisenberg term result in ionic displacements which cooperate with the 
purely electronic polarization, summing up to the appreciable value of 6 /iC/cm 2 . 

So far, we have discussed the prototypical case of HoMn03 ; however, as previously mentioned, 
the E-type is the magnetic ground state for many distorted manganites [97] and it is therefore 
interesting to investigate how the relevant properties (with a focus on polarization) change as 
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Table 1: Relevant calculated properties in HoMn03 . First two columns: Mn-O-Mn angles, 
broken down into values for the case of parallel (a p ) and antiparallel (a ap ) spin. Third to fifth 
columns: polarization values calculated when considering only the electronic polarization in the 
original centrosymmetric structure (P^), or only the PCM value upon structural relaxation 
(P^n M ) an d the total Berry-phase polarization for the relaxed ionic coordinates (P^f). Sixth 
to eighth columns: (3,3) components of the Born effective charge tensors, for Mn ions (Z*(Mn)) 
and the two inequivalent in-plane oxygens (Z*(O p ) and Z*(O ap )). 
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a function of the rare-earth |78j. Recall that the rare-earth cation has primarily the effect 
of increasing the octahedral GdFe03-like tilting as a result of reducing the ionic size when 
moving, say, from La to Lu; on the other hand, the Jahn- Teller-like distortions are weakly 
affected by the rare-earth atom |97[ [78] . The structural modifications (relative to the Mn-O- 
Mn angles) have in turn important consequences on the magnetic and dipolar order. As for 
the former, we have shown [78] that the first-nearest-neighbor ferromagnetic exchange-coupling 
constant progressively weakens upon decreasing the ionic radius, whereas the strong second- 
nearest-neighbor AFM exchange constant is more or less constant along the series. This implies 
the progressive change of the magnetic ground-state from A-type (in early rare-earth manganites) 
to E-type (in late rare-earth manganites), going through the intermediate region (R = Tb, Dy) 
where the spin-spiral occurs as ground state. What happens to polarization? To perform a 
complete investigation of the ferroelectric properties as a function of the octahedral tilting, we 
have imposed the E-type magnetic state on all the rare-earth manganites, irrespective of the 
actual magnetic ground-state. This is a typical example of a "computer-experiment": within 
DFT, at variance with real experimental samples, one can impose several different structural, 
electronic or magnetic configurations (not necessarily the ground states) to have clear insights 
on specific phenomena or to separate several competing effects. 

What we focus on here is the construction of Wannier functions (WF) [22\ [95] for the Mn e g , 
Mn t2 g and O p band manifolds and on the position of the WF center with respect to the 
relative ionic site. The difference between the polarization calculated according to the point- 
charge-model and via the Berry-phase approach is commonly referred to as the "anomalous" 
contribution to polarization. As such, it reflects somewhat the deviation from a purely ionic 
state or, equivalently, highlights the covalent character of the atomic bonds and, in turn, of the 
electronic structure. Moreover, we also recall that the polarization via the Berry-phase approach 
is equivalent to the sum of the displacement of the center of each WF from the position of the 
corresponding ion plus PCM contribution. The latter was shown [78] to be rather unaffected by 
the i?-ion, with a value P^ M ~ 2 fiC/cm 2 . 

In Fig. [7J we report the different contributions to the total polarization in the spin-up channel 
coming from the displacements of the WF centers for the Mn e g , Mn t2 g and O p, along with their 
sum (leading to the spin-up "anomalous contribution" ) . We note that Mn t^g states contribute 
in an opposite way with respect to Mn e g and O p states, the total P having the same sign 
as the two latter contributions. Moreover, it is quite clear that, whereas the O p and Mn t2 g 
depend relatively little on the rare-earth ions, the e g contribution is very sensitive to structural 
distortions. Indeed, for a hypothetical LaMnOs in the E-type spin configuration, there would be 
a total polarization (coming from twice the spin-up contribution shown in Fig. [7J plus the PCM 
term), summing up to a value greater than 10 //C/cm 2 ! This confirms the strong sensitivity 
of the e g states to the Mn-O-Mn angle: as reported in Ref. [78], the hopping integral strongly 
decreases when moving from La to Lu, consistent with a progressively reduced band width. 
Whereas promising ways to increase P would appear in the early rare-earth manganites (but 
where unfortunately the magnetic ground state is the (paraelectric) A-type AFM), the total 
polarization seems pretty much "saturated" to a value of the order of 6 /iC/cm 2 in going from 
Ho to Lu. 

We would like to comment now on the comparison with experiments. First of all, we remark 
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that several problems exist with the experimental synthesis of the late R manganites: indeed, 
the stable structure is hexagonal, not orthorhombic [98|, I97|. Modern growth techniques, such 
as high-pressure high-temperature synthesis, can do the job and synthesize ortho-manganites 
for late rare-earths, leading however not to single-crystals but rather to polycrystalline samples. 
This poses problems for the exact evaluation of ferroelectric polarization, due to possible different 
orientations of the polarization vector in the polycrystalline grains. To our knowledge, there 
exists several values in the literature. Lorenz et al. [99] reported P ~ 0.001 fiC/cm 2 for HoMn03, 
i.e. a value smaller by two or three orders of magnitudes than our ab-initio estimates. On the 
other hand, a much larger value was recently reported in AFM-E TmMnOs [100] : a lower bound 
of (unsaturated) polarization of about 0.15 //C/cm 2 was measured, in much better agreement 
with our theoretical values. This is especially so, since Pomjakushin et al. [100] suggested that 
the threshold of 1 //C/cm 2 could be easily achieved in the case of single crystals. In this respect, 
we would also like to remark that the values discussed so far are calculated within a bare DFT 
approach. It is however well known that DFT fails in accurately modelling strong correlation 
effects, which might occur in manganites. However, the inclusion of an Hubbard- like correction 
according to the so-called LSDA+C/ approach for Mn d states in HoMnC>3 , lead to values of 
the polarization all larger than 1-2/xC/cm 2 for U < 8 eV. Recently, in Ref. [101] . the authors 
reported a theoretical model in the context of electromagnon excitations in RM11O3. One of the 
outcome was the estimate of the polarization in E-type manganites based on optical absorption 
data measured for TbMnOs in the spiral-phase: P was found to be of the order of 1 //C/cm 2 , 
therefore large and compatible with our theory estimates. Though some controversy is still 
present, there are more and more confirmations that the polarization in E-type is much higher 
than in the spiral phases studied so far, consistently with the generally accepted argument that 
magnetostrictive effects in the symmetric Heisenberg-like exchange should be stronger than in 
the antisymmetric DM part. 

3.1.2 Half-doped manganites: Lao.sCao.sMnOa 

Hole-doped manganites {i.e. Ai- x B x MnO^ where A = La, Pr, . . . and B = Ca, Sr, . . . ) show 
a rich physics, with exciting phenomena ranging from charge-ordering to half-metallicity, from 
colossal magnetoresistance to exotic phase diagrams, from orbital-ordering to metal-insulator 
transitions. We will here discuss the possibility that hole-doped manganites, with a hole- 
concentration x ~ 0.5, might also become ferroelectric and, therefore, multiferroic. 

Lao.sCao.sMnOa (denoted in the following as LCMO) is a very complex system from many 
points of view (electronic, structural, magnetic, etc.) and, despite the many decades of work 
since the first seminal paper [102] . its properties have not been clearly elucidated. In particular, 
even the exact ionic coordinates and related symmetries are still debated. Two main models 
have been proposed so far: a) the first one, proposed by Radaelli et al. [103] . is based on 
a site-centered charge-ordered (SC-CO) Mn 3+ /Mn + checkerboard arrangement in the Mn02 
plane (see Fig. [8b), in which the octahedron around Mn 3+ is Jahn- Teller- like distorted, whereas 
the octahedron around Mn 4+ is rather regular; b) the second one, proposed by Rodriguez et 
al. [104] and referred to as a bond-centered charge-ordered (BC-CO), is based on a structural 
dimerization of Mn ions (all in a d A configuration). This leads to a peculiar OO: at variance 
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Figure 7: Different up-spin contributions to the "anomaious" term in the polarization (in 
//C/cm 2 ) as derived from WF centers: Mn e g (red), Mn t2 g (blue), O p (green) and total 
(magenta) as a function of the rare-earth ion (R = La, Nd, Sm, Ho, Lu). 
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Figure 8: a) Checker-board arrangement of Mn 3+ and Mn 4+ in the Mn02 plane in the SC-CO 
structure. The AFM-CE magnetic configuration is shown by double zigzag up (black arrows) 
and down (red arrows) spin chains, b) Sketch of the 9 rotation: the spins on two neighboring 
Mn atoms in the up-spin chain are rotated clockwise by 9 (green arrows), along with two 
corresponding spins on neighboring Mn in the down-spin chain rotated clockwise by 9 (blue 
arrows), c) The schematic orbital-ordering in the SC-CO structure: ideally, there should be 
an elongated Jahn- Teller-like e g orbital centered on the Mn 3+ site and no-e 9 -like charge on the 
Mn 4+ site, d) The schematic 00 in the BC-CO structure: the two Mn ions in the dimer show 
their e g orbitals oriented one towards each-other. ZP units {i.e. two Mn and the O in-between) 
are highlighted by ellipses. 

with the staggered 00 previously mentioned for LaMnOs, here the filled Mn e g orbitals in the 
dimer point one towards each other. With respect to the mother compound, LaMn03, there is 
one extra-hole every two Mn: the (spin-polarized) hole is believed to be located on the central 
O in between the two Mn. This peculiar unit (formed by two Mn and the O in between) is often 
referred to as "Zener-polaron" (ZP) [1051 1106] . after the Zener double exchange mechanism 
which should be enhanced here (see Fig. [8H) . 

As far as the magnetic spin-configuration is concerned, the so-called CE-type AFM (i.e. double 
zig-zag spin chains in the Mn02 plane, cfr. Fig. [Hh.) has been proposed as ground-state. 

We will here focus on two different mechanisms which might lead to improper ferroelectricity in 
LCMO: 

• The first one is based on breaking inversion symmetry in the spin-chains through a rotation 
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(by an angle 9) of the spins on two nearest-neighbor Mn in the up zigzag chain, along with 
a corresponding rotation of two spins in the down spin-chain (c/r. Fig. [8b), so as to keep 
a global AFM character. This follows the theoretical proposal put forward by Efremov et 
al. [107 \, who first suggested the possibility of multiferroicity in manganites. According 
to Ref. |107j . such rotation should progressively lead from a fully SC-CO (in the ideal 
CE-type, 9 = 0°) to a fully BC-CO for 9 = 90° (where the dimerization process driven by 
spin ordering is maximized). Efremov et al. predicted that, in both the extreme cases, 9 = 
0° and 9 = 90°, the polarization should vanish: for 9 = 0°, the checkerboard arrangement 
should be fully centrosymmetric (both structurally and electronically), whereas for 9 = 90° 
the Mn should not show any charge-disproportionation. However, for in-between values of 
9, the intermediate SC-CO/BC-CO should lead to a small charge-disproportionation and, 
therefore, to inequivalent Mn (at variance with the ZP state and reminiscent of the site- 
centered CE-type). In this case, inversion symmetry would be broken by spin-dimerization, 
therefore paving the way to ferroelectricity; 

• The second mechanism occurs in the structure experimentally proposed by Rodriguez et 
al. [104J. The related unit cell shows a "structural' Mn-Mn dimerization and implies 
a realization of a BC-CO, not invoking (non-collinear) magnetic mechanisms as in the 
previous case, but rather thanks to electronic rearrangement — such as OO — following 
the structural distortions. Still, in this case, our mechanism for multiferroicity is once 
more a (collinear) magnetically induced mechanism based on the inequivalency of some 
specific oxygen atoms, as will be detailed below. 

Due to the large unit-cell (80 atoms, needed to simulate the CE-type AFM ordering, along with 
a checkerboard arrangement of La and Ca cations) and the need of non-collinear spin magnetism 
(needed to simulate finite values of 9), the computational cost of these simulations is very high. 
For this reason, the ionic positions were not optimized within DFT, but were rather kept frozen 
in the structure proposed either by Radaelli [103] or by Rodriguez [104] . labelled in what follows 
as LT-M or by LT-O, respectively. Unfortunately, the lack of ionic minimization forbids any 
DFT prediction of the actual structural and magnetic ground-state from total-energy arguments; 
this calls for future studies. From our calculated values for unrelaxed structures, it seems that 
the CE with SC-CO is the phase showing lowest total energy; however, for example, the SC-CO 
state with a rotation 9 = 45°, is higher in energy by only ~ 4 meV/Mn. One can therefore 
conjecture that, in real samples, there might be a coexistence of nanoscale regions with different 
magnetic structures (i.e. with zero and finite 9 values). 

Before discussing the relevant ferroelectric properties, let us mention some general features 
in the electronic structures of the LT-M and LT-0 systems, both in the CE-type AFM spin 
configuration (i.e. 9 = 0). In Fig. [9] we show the isolines of the electronic charge plotted in the 
energy region where the Mn e g states are located. It is clear that in the LT-M (Fig. [9^) the 
shape of the e g electronic cloud, centered on the "nominal" Mn 3+ , is markedly elongated towards 
the neighboring Mn 4+ with parallel spins. On the other hand, the Mn 4+ show a very isotropic 
distribution of the charge. The situation is different in the LT-0 structure (Fig. [9b), where the 
OO clearly shows the e g orbitals forming "dimers" with their charge distribution pointing one 
towards the other, as driven by the underlying ionic configuration. Let us mention a note on the 
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Figure 9: Isolines of the e g charge in a) the LT-M SC-CO and b) the LT-0 BC-CO structures. 
Red (blue) lines marks the minimum (maximum) charge, through the intermediate green lines. 
In a), black arrows mark the Mn spin directions. In b), ZP are highlighted. 

CO: consistently with previous reports, the actual charge-disproportionation in LCMO within 
DFT is of the order of only 0.1-0.2 electrons in the LT-M SC-CO, at variance with the ideal 
situation of "full" charge disproportionation, where the e g electron cloud should be completely 
distributed around the Mn 3+ , with no-charge on the Mn 4+ . In this sense, the calculated 00 
in the LT-M (c/r. Fig. [9b.) is different from the nominal situation (c/r. Fig. [8b) with clear 
signatures of e g charge also around the Mn 4+ . We remark, however, that the small charge- 
disproportionation detected in the LT-M structure becomes really negligible (<0.02 electrons) 
in the LT-0 BC-CO; this suggests that it is still meaningful to consider the LT-M — » LT-0 
transition as a corresponding SC-CO — > BC-CO transition. 

Let's now consider what happens in the LT-M structure upon increasing 8 from the initial 
zero- value: our calculated electronic structures (not shown, see Ref. |85j ) indicate a decreasing 
e g band-width and a related increasing band-gap. This can be rationalized by comparing the 
spin-arrangement with finite 6 with the original CE-type. Upon spin rotation, the e g electron 
- which could hop equivalently on the two nearest-neighbor Mn 4+ on both sides along the 
spin-chain in the CE-AFM phase — will now preferentially hop on the Mn 4+ which shows a 
parallel spin, since hopping in the other direction is prevented by the spin misalignment. This 
effect would rather lead to a decreased hopping and to a reduced e g band-width, at variance 
with our findings. However, one needs to consider that, for 6^0, there will be an increasing 
probability to hop on the neighboring spin-chain (prevented by opposite spin configuration in 
the CE-AFM phase). Overall, there will be therefore an increased hopping integral, consistently 
with our findings and with Hamiltonian- modelling studies, as well [85]. 
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For our purposes, the most important finding is that a finite 9 in the LT-M induces a rather 
large polarization, as shown in Table [2] (first line), with an increasing parabolic trend of P 
vs. 9. We note that this is a purely electronic polarization, since the ions are fixed in their 
centrosymmetric arrangement [103] , The Heisenberg-like symmetry breaking — as driven by 
spin-rotation — is therefore confirmed as an efficient tool to induce large ferroelectricity (recall 
that spin-orbit coupling and the related DM interaction is neglected in the present context). 

We will now focus on the calculated values of P in the LT-0 structure (see Table [2]) and start 
the discussion for the 9 = case. Without any magnetic ordering imposed and as determined 
experimentally, the LT-0 structure shows a P1\nm space group: this implies that some in-plane 
oxygens are structurally equivalent (shown in the same color in Fig. [10]) . due to the 1\ screw 
symmetry. However, when imposing the AFM-CE spin-configuration, the Oxygen equivalency 
is lifted: there is an alternation of O p bonded to two parallel spins and of O ap bonded to two 
antiparallel spins. This is sufficient to give rise to ferroelectricity in the direction shown in Fig. 1 101 
Remarkably, the induced polarization reaches the surprisingly large value of several /iC/cm 2 . 
To further verify that a magnetically-induced mechanism is the source of the ferroelectricity, we 
have also performed a 9-\ike rotation of the spin dimers, similar to the previous case of the LT-M 
(c/r. Fig. 0b)). In this case, upon spin-rotation, the inequivalency of the Oxygens crossing the 
2i axis is reduced. In the extreme situation, 9 = 90°, all the O atoms are now bonded to two 
Mn with perpendicular spins: in this configuration, they all look equivalent and the source of 
polarization vanishes. Indeed, DFT calculations confirm that this is the case (c/r. Table [2]). In 
summary, our DFT results (both from H0M11O3 and LCMO) offer a confirmation that the O 
inequivalency is an efficient handle to achieve and/or tune a large ferroelectric response. 

3.2 Problems and perspectives in Improper multiferroics 

As indicated by the huge interest in the last few years, magnetically-driven ferroelectrics, with 
ortho-TbMn03 taken as prototype, are with no doubt an exciting class of materials. However, 
there are a few bottlenecks which prevent their use in large-scale applications: i) their polariza- 
tion is generally very small (< 0.1 //C/cm 2 ); ii) their ordering temperature is very low (of the 
order of few tens of K); Hi) being globally antiferromagnets, their net magnetization is always 
zero (a ferromagnetic spin ordering alone cannot break inversion symmetry!). In this respect, 
we will certainly see some activity in future years to get rid of these problems. 

As shown in this review, at least point i) can be beautifully overcome when considering Heisenberg- 
like exchange-striction, as shown in E-type manganites. The ordering temperature of the latter 
is, however, extremely low (Tjy(HoMn03 ) ~ 26 K). One possibility to increase the ordering 
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Table 2: Berry-phase polarization (in //C/cm 2 ) calculated in the LT-M and LT-0 structure as 
a function of spin-rotation angles 9 (first line). 
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Figure 10: Atomic configuration in the Mn02 plane of the LT-0 structure: symmetry-equivalent 
atoms are marked in the same color. Note that red and black spheres mark Mn atoms: despite 
being symmetry-inequivalent, the two kinds of Mn are only marginally different from the elec- 
tronic point of view, with small differences in the Mn-0 bond- lengths (see Ref. [104] . ) . Hori- 
zontal lines mark the two 2i screw axes in the unit cell, crossing the O atoms (marked as grey 
and yellow). The blue arrows on the Mn ions denote the spin directions in the AFM-CE spin 
configuration: when considering the spin-directions, the grey atoms (structurally equivalent by 
symmetry) become electronically different: they are alternatively bonded to two parallel Mn 
spins and to two anti-parallel spins (see labels on two selected oxygens). 

temperature without losing the non-centrosymmetric Heisenb erg-like exchange-striction is to 
consider rare-earth nickelatcs [86] (for example, T/v(HoNi03) = 145 K, Tjy(LuNi03) = 130 K, 
etc.). Nickelates are rather complex materials, with several important issues still under debate, 
including the origin of their metal-insulator transition as well as their spin configuration. As 
for the latter, both non-collinear and collinear spin-arrangements have been put forward from 
neutron diffraction studies [108^ 1109] , In addition, nickelates show a charge-disproportionation: 
Ni ions, in the nominal 3+ valence-state, split into two groups of Ni 2+ and Ni 4+ [110] . This adds 
one degree of freedom to achieve ferroelectricity. For example, as suggested in Ref. [Ill] , one 
of the proposed magnetic configurations shows, along the [111] direction, a sequence of Ni 2+ - 
Ni 4+ -Ni 2+ -Ni 4+ as for charge-ordering and a sequence of till planes as for spin-ordering. The 
combination of spin and charge-ordering would break centrosymmetry, leading to a polarization 
along the [111] direction. Another spin-configuration, proposed by experiments, seems to be 
very similar to the E-type in HoMn03 , the only difference being the stacking of TMO2 (TM 
= Mn, Ni) planes: whereas the out-of-plane coupling is always AFM in HoMn03 , in nickelates 
there are MO2 alternatively coupled ferromagnetically and antiferromagnetically. However, the 
different out-of-plane stacking does not destroy the mechanism for polarization, induced in a 
way very similar to HoMn03 . Our preliminary calculations [86] show that the two mentioned 
collinear magnetic ground-states in monoclinic i?Ni03 (R = Ho, Lu) are basically degenerate 
(i.e. the differences in total energies are below our numerical uncertainty). Consistently with a 
Heisenberg-driven mechanism, both spin-configurations give rise to a large polarization (of the 
order of few /uC/cm 2 ) along different directions, suggesting nickelates as a new and interesting 
class of magnetically-driven multiferroics. 
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Going back to the bottlenecks mentioned above, point iii) might be overcome by considering 
magnetite. In this review, we have discussed so far a few examples where spin-ordering is 
a necessary ingredient to break inversion symmetry. However, there are materials in which 
the polarization is induced purely by charge-ordering, such as LuFe204 and Fes04 below the 
Verwey transition temperature (i.e. corresponding to the metal-insulator transition, Ty ~ 
120 K). In magnetite, the spin-arrangement is ferrimagnetic (i.e. tetrahedral and octahedral 
Fe sites show up and down spin, respectively). The role of magnetism, however, does not 
seem to be relevant for polarization. Magnetite is a complex and controversial system: the 
Fe 2+ /Fe 3+ charge ordering pattern on octahedral iron sites is still under debate [1121 1113] . 
However, the Cc symmetry has been proposed by diffraction studies and confirmed from first- 
principles to be the ground state [114] . In the Cc case, octahedral Fe sites, form a corner-sharing 
tetrahedron network: 75% of the tetrahedra show the so-called "3:1" pattern (meaning that, 
in each tetrahedron, 3 sites are Fe 2+ and one is Fe 3+ or vice versa), whereas 25% show a 2:2 
pattern (meaning that 2 sites are Fe 2+ and two are Fe 3+ in the tetrahedron). It happens that the 
Cc is non-centrosymmetric; indeed, our DFT calculations [HI] show the polarization induced by 
charge-ordering to be of the order of few //C/cm 2 , suggesting magnetite to be the first improper 
multiferroic known to mankind. 

4 Summary and Conclusions 

In summary, we have presented some examples which show the power of DFT-based methods in 
the field of multiferroic materials. This includes: i) rationalizing experimental observations in 
known multiferroics, ii) designing new (artificial) multiferroics with optimized properties (larger 
ferroelectric polarization, strong ferromagnetism, higher ordering temperatures, etc.), and iii) 
proposing and quantifying novel microscopic mechanisms, based on electronic degrees of freedom, 
which potentially lead to ferroelectricity in magnetic transition metal oxides. 

It is apparent that the field of proper magnetic ferroelectrics has a relatively long history: many 
of these materials have already been studied in the 1960s or later, but have only recently been 
rediscovered. Due to substantial advancements in experimental synthesis and characterization 
techniques on one side, and the availability of powerful computational methods together with new 
theoretical approaches on the other side, substantial progress in understanding these materials 
has been achieved during recent years. Similar to the the case of non-magnetic ferroelectrics, 
first-principles calculations have shown a remarkably high degree of accuracy, reliability, and 
predictive capability for the class of proper multiferroics. Nevertheless, many open questions 
still remain, in particular how to achieve large polarization, large magnetization, and strong 
magneto-electric coupling above room temperature, or what mechanisms for coupling between 
magnetic and ferroelectric properties do exist in these materials. 

On the other hand, the field of DFT calculations for improper multiferroics is only a couple of 
years old. As such, it is not clear at the moment how accurate the predictive capabilities of cur- 
rent DFT approaches are for relevant quantities such as structural or electronic properties and, 
most importantly, polarization. On the experimental side, the synthesis of some compounds (i.e. 
as shown for ortho-manganites with late rare-earth ions) is not under full control, making the 
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theory-experiment comparison rather complicated. On the modelling side, the role of electronic 
correlations (where DFT often shows its limits) is certainly more relevant in improper than in 
proper magnetic ferroelectrics. In this respect, future developments on the theory side (i.e. 
invoking novel exchange-correlation functionals to better describe many-body effects) are desir- 
able. As such, a strong interaction with the experimental and model-Hamiltonian communities 
active in the field, as well as the extension of DFT studies to a much larger set of materials 
(showing different microscopic mechanisms or simply different chemical, structural, or electronic 
properties), will be necessary to achieve a satisfactory qualitative and quantitative description 
of the complex physics at play in improper multiferroics. 
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